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Abstract 

This paper describes the calibration procedure for the drift tubes of the CMS barrel 
muon system and reports the main results obtained with data collected during a high 
statistics cosmic ray data-taking period. The main goal of the calibration is to deter- 
mine, for each drift cell, the minimum time delay for signals relative to the trigger, 
accounting for the drift velocity within the cell. The accuracy of the calibration pro- 
cedure is influenced by the random arrival time of the cosmic muons relative to the 
LHC clock cycle. A more refined analysis of the drift velocity was performed during 
the offline reconstruction phase, which takes into account this feature of cosmic ray 
events. 
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1 Introduction 



The Compact Muon Solenoid (CMS) d is a general-purpose detector whose main goal is to 
explore physics at the TeV scale, by exploiting the proton-proton collisions provided by the 
Large Hadron Collider (LHC) at CERN. 

CMS uses a right-handed coordinate system, with the origin at the nominal collision point, 
the x-axis pointing to the center of the LHC, the y-axis pointing up (perpendicular to the LHC 
plane), and the z-axis along the anticlockwise-beam direction. The polar angle, 9, is measured 
from the positive z-axis and the azimuthal angle, <p, is measured in the x-y plane. 

The central feature of the Compact Muon Solenoid apparatus is a superconducting solenoid, of 
6 m internal diameter, providing a field of 3.8 T. Within the field volume are the silicon pixel 
and strip tracker, the crystal electromagnetic calorimeter (ECAL) and the brass/scintillator 
hadron calorimeter (HCAL). Muons are measured in gas-ionization detectors embedded in the 
steel return yoke. In addition to the barrel and endcap detectors, CMS has extensive forward 
calorimetry. 

The barrel muon system [3J is divided in five wheels. Every wheel is composed of 12 sectors, 
each covering 30° in azimuth, as shown in Fig.JT] 



Figure 1: Schematic representation, in the x — y plane, of the chamber positions within a wheel 
of the muon barrel system of the CMS experiment. The labels and the numbers of the muon 
stations are shown. Because of mechanical requirements, the top and bottom MB4 sectors are 
split in two distinct chambers. 

Each sector contains four stations equipped with Resistive Plate Chambers (RPC) and Drift- 
Tubes (DT) chambers. The four DT chambers are labeled MB1, MB2, MB3, and MB4 going 
inside-out. In total there are 5 wheels * (3 stations * 12 sectors + 1 station * 14 sectors) = 250 
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DT chambers. The chambers are interleaved with the steel return yoke of the magnet and are 
composed of three groups, called "super-layers" (SL), of four staggered layers of independent 
drift cells, for a total of about 172 000 channels. A schematic representation of a chamber is 
shown in Fig. |2](left). 

The chamber volume is filled with a Ar(85%)/C02(15%) gas mixture, kept at atmospheric 
pressure. Two of the super-layers have the wires parallel to the beam direction and measure 
the rep coordinate, the other super-layer has wires perpendicular to the beam direction and 
measures the z coordinate. The chamber provides a measurement of a track segment in space. 
The outermost station is equipped with chambers containing only the two r<p super-layers. The 
basic element of the DT detector is the drift cell, illustrated in Fig. [2] (right), where the drift lines 
and isochrones are represented. All chambers were operational, fully commissioned, and the 
number of problematic channels were less than 1 %. 

The DT system is designed to provide muon track reconstruction, with the correct charge as- 
signment up to TeV energies, and first-level trigger selection. It also yields a fast muon iden- 
tification and an accurate online transverse momentum measurement, in addition to single 
bunch-crossing identification with good time resolution. The good mechanical precision of the 
chambers allows the track segments to be reconstructed with a resolution better than 250 f/m 
0. 

A fundamental ingredient of the DT system is the calibration, which is used as input to the 
local hit reconstruction, within the drift cells, and thereby influences the precision of the track 
reconstruction. This paper describes in detail the DT calibration procedures, to obtain reliable 
calibration constants, and presents results obtained from the extensive commissioning run with 
cosmic ray events performed in Autumn 2008, the Cosmic Run At Four Tesla (CRAFT) [4], in 
preparation for LHC running. 




Figure 2: Left: Schematic view of a DT chamber. Right: Section of a drift tube cell showing drift 
lines and isochrones. The voltages applied are +3600 V for wires, +1800 V for electrode strips, 
and — 1200 V for cathode strips. 

The data sample used for the calibration process and the trigger conditions are summarized 
in Section 2. The main characteristics of the DT calibration process are described in Section 3. 
The process consists in the determination of the inter-channel synchronization, described in 
Section 4, the analysis of noisy channels, treated in Section 5, and the calculation of the time 
pedestals and the drift velocity, described in Sections 6 and 7, respectively. The DT calibration 
workflow, including the monitoring of the conditions, is performed within the CMS computing 
framework, as described in Section 8. Finally, a more refined analysis of the drift velocity in the 
muon system, within the offline reconstruction process, is presented in Section 9. 
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2 Cosmic Ray Event Trigger and the Data Sample 

The long cosmic ray data taking in 2008 with and without magnetic field allowed a detailed 
study of the DT drift properties and an improved understanding of the calibration constants. 
About 270 million events were collected with a 3.8 T field inside the solenoid magnet. In this 
configuration, the radial component of the magnetic field in the DT chamber positions does not 
exceed 0.8 T. 

The DT system was the primary trigger source for most of the collected events. The local trigger 
[3j was designed to operate with collisions taking place at the center of the CMS detector, 
and it is performed searching the ^-matching of hits in each chamber. This is achieved using 
dedicated hardware, which configures the expected track paths from one chamber to another. 

Due to the different origin, direction, and timing of the cosmic rays, as compared to muons 
from proton-proton collisions, dedicated adjustments were needed to properly configure the 
DT trigger for high efficiency during CRAFT. This required relaxing the extrapolation algo- 
rithm with a particular configuration of the DTTF (DT Track Finder) as explained in more detail 
in Ref. |5|. Therefore, during data-taking with cosmic rays, the LI trigger was generated by the 
coincidence in time of two segments in two stations of the same sector, or adjacent sectors, and 
a rate of about 240 Hz was provided to the Global Muon Trigger. 

About 20 million events, out of the 270 million collected during CRAFT, were used for the 
calculation of the calibration constants. They have been chosen from stable runs where most 
of the DT system was operational. No quality cuts are, in principle, necessary to perform the 
calibration. However, in order to have a clean sample of muons, a transverse momentum cut 
of 7 GeV was applied. 

3 The Calibration Process 

Charged particles crossing a DT cell produce ionization electrons in the gas volume. The de- 
termination of the relationship between the arrival time of the ionization signal and its spatial 
deposition is the primary goal of the calibration task, which leads to the extraction of the drift 
times and drift velocities. 

The arrival time of the ionization signal is measured using a high performance Time to Digital 
Converter (TDC) [6|. This is the main building block of the read-out boards of the DT system. 
It is a multi-hit device in which all hits within a programmable time window, large enough to 
accommodate the cell maximum drift time, are assigned to each Level 1 Accept trigger. The 
drift time is directly obtained from the time measured by the TDC, after subtracting a time 
pedestal which contains contributions from the latency of the trigger and the propagation time 
of the signal, within the detector and the data acquisition chain. The first goal of the calibration 
procedure is, therefore, to determine the time pedestals, as described in Section 6. The expected 
precision of the time pedestal calibration during the cosmic ray data-taking is limited by the 
arrival time distribution of cosmic rays which is flat within the clock cycle and it is of the order 
of 25 ns/ a/12. 

The other relevant quantity for the DT calibration process is the effective drift velocity. It de- 
pends on many parameters, including the gas purity and the electrostatic configuration of the 
cell, the presence of a magnetic field within the chamber volume, and the inclination of the 
track. The parameters connected to the working conditions of the chambers are monitored 
continuously [3j: the high voltage supplies have a built-in monitor for each channel; the gas 
is at room temperature and its temperature is measured on each preamplifier board inside the 
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chamber; the gas pressure is regulated and measured at the gas distribution rack on each wheel, 
and is monitored by four further sensors placed at the inlet and outlet of each chamber. The 
adequacy of the flow sharing from a single gas distribution rack to 50 chambers is monitored 
at the inlet and outlet line of each individual chamber. A possible leakage in the gas line can be 
sensed via the flow and/ or the pressure measurements. 

Five small gas chambers, one for each wheel, are used to measure the drift velocity in a volume 
of very homogeneous electric field, located in the accessible gas room adjacent to the cavern, 
outside of the CMS magnetic field. Each of these chambers, called Velocity Drift Chambers 
(VDC) EL is able to selectively measure the gas being sent to, and returned from, each in- 
dividual chamber of the wheel thus providing rapid feedback on any changes due to the gas 
mixture or contamination. During the CRAFT data-taking period only one such chamber was 
used. 

No noticeable variation of the parameters described above is expected among different regions 
of the spectrometer. However, the magnetic field and the track impact angle may vary substan- 
tially from chamber to chamber, as they occupy different positions in the return yoke. 

Two methods for calculating the electron drift velocity in a DT cell are presented here. They 
both assume a constant drift velocity in a given chamber. The first method, discussed in Sec- 
tion 7, is based on the measurement of the effective drift velocity using the mean-time tech- 
nique, which computes the velocity value at the super-layer granularity level. The second 
method, discussed in Section 9, relies on the muon track fit, which determines track-by-track 
the time of passage of the muon and the drift velocity as additional free parameters of the fit, 
together with the track position and inclination angle. The assumption of a constant drift ve- 
locity is considered a good approximation because the magnetic field in the chamber volume 
is usually low and approximately homogeneous. A third method based on a parameterization 
of the drift velocity as a function of the drift time, the magnetic field, and the muon trajectory 
is discussed in Ref. 13. 

The detailed drift velocity analysis, described in Section 9, also reveals non-linear effects in 
the innermost stations (MB1 chambers) of the barrel external wheels (Wheel +2 and Wheel -2), 
where the strongest radial magnetic field component, of about 0.8 T, is present. 

The DT calibration process also depends on the different signal path lengths to the read-out 
electronics (called inter-channel synchronization time) and on the list of noisy channels, as will 
be described in the following sections. 

4 The Inter-Channel Synchronization 

The inter-channel synchronization is calculated for each read-out channel of each chamber, in 
order to correct for the different signal path lengths of trigger and read-out electronics. This 
is a fixed offset, since it only depends on cable /fiber lengths, and it does not need to be re- 
calibrated very often. Nevertheless, it is useful to frequently redo its calibration, to monitor the 
correct behavior of the front-end electronics. The inter-channel synchronization is determined 
by test-pulse calibration runs. The design of the data acquisition system allows such runs to 
be taken during the normal physics data-taking, by exploiting the collision-free interval of the 
LHC beam structure, called "abort gap". 

During special calibration runs, a test-pulse is simultaneously injected in four channels of a 
front-end board, each one from a different layer of a super-layer, simulating a muon crossing 
the super-layer. To perform the scanning of the entire DT system in only 16 cycles, the same 
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test-pulse signal is also distributed to other four-channel groups, 16 channels apart. 

The so-called fo calibration consists in determining, for each DT channel, the mean time and 
the standard deviation of the test-pulse. In the calibration procedure, the events are split in 
two samples: the first is used to compute the average value, within a full chamber, of the signal 
propagation time from the test-pulse injector to the read-out electronics; the second is used to 
calculate, for each individual channel, the difference between the time of its test-pulse signal 
and the average value of the chamber. 
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Figure 3: Inter-channel synchronization constants calculated from a test-pulse run. The results 
are shown for three representative layers, belonging to each of the three super-layers of cham- 
ber MB3 in Sector 9. The step-function shape reflects the grouping of channels among different 
front-end boards. 

Figure [3] shows an example of a distribution of fo constants for representative layers of the 
three super-layers of a chamber, as a function of the channel number. The other layers show 
very similar to values. The fo synchronization correction is always below 10 ns (1 TDC count 
corresponds to 0.78 ns). The standard deviation is about 1 ns, for all channels. This is com- 
patible with the precision of the electronic chain. These corrections correspond to the distance 
between the front-end boards, located inside the chamber volume, and the read-out boards. 
Wires connected to a given front-end board belong to cells adjacent to each other, in the super- 
layer, and have approximately the same distance up to the read-out boards. This leads to the 
step-function shape seen in Fig. [3], more pronounced in the super-layers 1 and 3. 

5 Noise Analysis in the DT Chambers 

On the basis of systematic studies performed during several commissioning phases of the DT 
detector, a cell is defined as noisy if its hit rate at operating voltage, counting signals higher 
than a common discriminator threshold of 30 mV, is higher than 500 Hz. 
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Figure 4: Distribution of the cell noise rate for different data-taking conditions: with and with- 
out magnetic field; with and without Cathode Strip Chambers (CSC), Resistive Plate Chambers 
(RFC), and Pixels (PIX). 



The number and geometrical distribution of noisy DT channels have been studied, in particular, 
during a commissioning period without magnetic field and having the detector wheels sepa- 
rated from each other. In this section we describe the results of the noise analysis performed 
using the cosmic ray events collected during CRAFT. With respect to runs using random trig- 
gers, the noise analysis based on normal data taking runs has the advantage of reflecting the 
detector operation in more realistic conditions. 

The first aim of the noise studies is to check the stability of the number of noisy cells in different 
conditions of the CMS detector. For all runs analyzed, the number of noisy cells is around 
0.01 % of all DT channels. The rate of noise hits per cell is shown in Fig. |4j for a number of 
representative runs, with and without magnetic field, and with different sub-detectors included 
in the acquisition. The number of cells with a hit rate higher than 500 Hz is very small. Detailed 
information on the noise rate observed for representative runs, with different configurations of 
CMS detectors included in the data acquisition, is shown in Table [T] 

An average noise rate of ~ 4 Hz is observed in the DT system, essentially insensitive to the 
magnetic field and to the status of nearby sub-detectors. In addition, it has been observed that 
around 50 % of the noisy cells remain noisy for long data-taking periods. 

Studies have also been made concerning the position dependence of the noisy cells, within 
wheels and chambers. As seen on Table [TJ most of the noisy channels are located in the in- 
nermost chambers (MB1), where the internal cabling is more complex, because of the reduced 
space. In Fig. U the noisy cell distribution is shown as a function of the wire number. The 
noisy cells appear concentrated in the regions of the super-layer close to the wire boundaries, 
which are different depending on the number of wires present in each chamber type (about 50 
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Table 1: Number of Noisy Cells in each chamber type (MB1, MB2, MB3, MB4) and average 
noise rate for some representative runs of different data-taking conditions. The results from the 
CRUZET (Cosmic RUn at ZEro Tesla) commissioning period are also shown, for comparison. 
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for MB1, 60 for MB2, 70 for MB3, and 90 for MB4). The peaks also reflect the position of the 
connectors distributing the HV inside the super-layer, which generate some electronic noise in 
their proximity. 
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Figure 5: Number of noisy channels as a function of the wire number, for the data-taking peri- 
ods mentioned in Table [l] Each distribution corresponds to a different run and shows the noisy 
cells observed for chamber types MB1, MB2, MB3 and MB4. 



The observed fraction of noisy cells (0.01 %) and the average noise rate (~ 4 Hz) in the full DT 
system are too low to affect the digitization efficiency or the trigger rate. It is important, how- 
ever, to exclude the noisy cells in the calibration process described in the following sections. 
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6 The Time Pedestal Calibration 

6.1 Computation of the Calibration Constants 

The time pedestal calibration is the process which allows the extraction of the drift time from 
the TDC measurement. For an ideal drift cell, the time distribution coming from the TDC (£tdc) 
would coincide with the distribution of drift time Adrift)/ an d would have a box shape starting 
from a null drift time, for tracks passing near the anode, up to about 380 ns, for tracks passing 
near the cathode. 

Experimentally, some non-linear effects related to the electric field distribution inside the drift 
cell have to be considered in the response of these cells; they are enhanced by the track incli- 
nation and by the presence of the magnetic field. In addition, different time delays, related 
to trigger latency, and different cable lengths of the read-out electronics, also contribute to the 
TDC measurements. The time measured by the TDC, £tdo can be expressed as 

^TDC = to + ^TOF + ^prop + ^Ll + ^drift / (1) 

where 

• fo is the inter-channel synchronization used to equalize the response of all the chan- 
nels at the level of each chamber, as described in Section 4; 

• ^tof is the Time-Of -Flight (TOF) of the muon, from the interaction point to the cell, 
in the case of collision events. In the case of cosmic events, this quantity cannot be 
defined because the time pedestal is an average of the arrival time of cosmic muons 
relative to the clock cycle; 

• £prop is the propagation time of the signal along the anode wire; 

• £li is the latency of the Level-1 trigger; 

• tdrift is the drift time of the electrons from the ionization cluster to the anode wire 
within the cell. 

The main goal of the calibration is the calculation of the time pedestal, £ t rig/ which is dominated 
by the time delay caused by the LI trigger latency: 

ftrig = f TOF + £prop + ^Ll • (2) 

The value of ftrig is extracted for each super-layer directly from the fjDC distribution, referred 
to as Time Box, after subtracting the noisy channels and correcting for the inter-channel syn- 
chronization. Figure [6] shows a Time Box measured during a CRAFT run, for one super-layer. 

The value of £ t rig is the turn-on point of the Time Box distribution. It is computed by fitting 
the rising edge of the distribution to the integral of a Gaussian function, as illustrated by the 
continuous line in Fig. [6j The procedure is applied at the super-layer level and is described in 
more detail in Ref. |T0|. 

The main quantities calculated by the fit are the inflexion point of the rising edge, T mean , and 
its standard deviation, T S i gma , which represents the resolution of the measurement. Figure [7] 
shows the distributions of T mea n and T S i gma measured, in a CRAFT run with B = 3.8 T, for 
the innermost rep super-layer of a representative wheel, as a function of chamber type and 
sector. Similar results were obtained for the other wheels. Approximately constant values are 
observed for chambers of the same type and for all the wheels. The periodic structure seen in 
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Figure 6: Distribution of the signal arrival times, recorded by the TDC, for all the cells of a single 
super-layer in a chamber, after the cell-to-cell equalization based on the test-pulse calibration. 
The continuous line indicates the fit of the Time Box rising edge to the integral of a Gaussian 
function. 

the T mean distribution, Fig.[7|(top), reflects the time-of-flight of the cosmic muon from the upper 
sector to the lower sector. Indeed, the events contributing to the calculation of T mean and T S i gma 
can be triggered by the upper or lower sectors. The events triggered by the top (bottom) sectors 
may also be detected by other non-triggering sectors, having a less precise time pedestal and, 
consequently, leading to a less precise determination of these quantities. Different runs during 
the entire CRAFT period have been analyzed, and a stable performance of the whole DT system 
has been observed. As expected, no dependence on the magnetic field strength was observed. 

The time resolution distribution, Fig. [7] (bottom), indicates the precision which the calibration 
procedure can reach with cosmic rays. A standard deviation of ~ 10 ns is observed for all 
super-layers in all wheels, except in the vertical sectors, where the number of events is limited 
and the muon crossing angles are large. The time resolution precision is limited mainly by the 
random arrival time of cosmic muons relative to the clock cycle. Furthermore, the resolution 
in Sector 1 is systematically worse than in Sector 7 because the trigger cables that distribute the 
Level 1 accept signal are longer and, therefore, generate larger skews in the signal transmission. 

After the determination of T mean and T S i gma , the time pedestal, ftrig/ is estimated as 



^trig — T mean k • T s jg ma . (3) 

The k factor is evaluated by minimizing the position residuals, using the local reconstruction 
of track segments within chambers. After a few iterations, a k factor of 0.7 was computed 
for the CRAFT data and was applied to all super-layers. The position residuals were then 
recalculated and a final correction to the time pedestals was computed dividing the remaining 
offsets observed in the residual distributions by a constant drift velocity (54.3 /mi/ns). The 
final itrig constants were stored in a database, as described in Section 8. 
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Figure 7: Mean (top) and standard deviation (bottom) of the fitted inflexion point of the Time 
Box rising edge, for the innermost r<p super-layer for a representative wheel. The triggering 
sectors (3, 4, 5 and 9, 10, 11) are synchronized among each other. The sectors with vertical 
chambers (sectors 1 and 7) detect much less cosmic ray muons, leading to a poorly defined 
rising edge and a less accurate calibration. 
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6.2 Validation of the Calibration Constants 

Once the f trig constants are computed, the calibration process proceeds with the validation step, 
which consists in studying the effect of these constants on the reconstruction algorithm. The 
analyzed quantities are the residuals computed, layer by layer, as the distance between the 
hit and the intersection of the 3D segment with the layer plane. A complete description of 
the local reconstruction procedure is given in Ref. IfTTl . To correct for the propagation time 
along the wire, the reconstruction of the segment is done in a multistep procedure. First the 
reconstruction is performed in the rty and z projections independently. Once two projections 
are paired and the position of the segment inside the chamber is approximately known, the 
drift time is corrected for the propagation time along the wire and for the TOF within the 
super-layer, and the 3D position is updated. 

The mean values of the residual distributions calculated for the innermost rep super-layer for 
a representative wheel are shown in Fig. [E] (top). A systematic offset with respect to the ori- 
gin is observed compatible with the systematic delay between the arrival time of the cosmic 
muon events and the clock cycle. The standard deviations of the fit to the residuals, shown 
in Fig. [8] (bottom) and in the range 400-600 }im, represent the spatial resolution obtained with 
the calibration process, a factor of two worse than the nominal resolution of about 250 jim (3|. 
The difference is caused by the spread of the muon arrival times inside the 25 ns time window 
associated with the LI trigger. This dilution will not occur with LHC collision data. 



7 The Drift Velocity Calibration 

The aim of the drift velocity calibration is to find the best effective drift velocity in each region 
of the DT system. In order to be consistent with the f trig calculation, described in Section 6, the 
drift-velocity calibration is computed with a super-layer granularity. 

The calibration algorithm is based on the mean-time technique described in detail in Ref. IfTOl - 
In this method, the maximum drift time in a cell, T max , is calculated considering nearby cells in 
three adjacent layers and using a linear approximation to determine the average drift velocity. 
As an example, Fig. [9] shows the simplest pattern of a muon crossing a semi-column of cells, 
together with the equations used to calculate T max . In general, T max depends on the track incli- 
nation and on the pattern of cells crossed by the track. Taking into account these dependencies, 
a spread of about 28 ns has been observed in the calculation of T max from CRAFT data. 



The effective drift velocity can be estimated assuming a linear space-time relationship, 

eff ^semi-cell 

"drift ~~ ^ t \ ' 
^- 'max 

where L semi . ce u = 2.1 cm is half the width of a drift cell. 



(4) 



Drift velocities measured for each chamber/ sector and for a representative wheel (other wheels 



give similar values) are shown in Fig. 10 for two CRAFT runs, one without (top) and one with 
(bottom) magnetic field. The drift velocity has approximately a constant value of 54.3 jinx/ ns, 
although with some systematic deviations, caused by limitations of the calibration procedure 
applied to cosmic ray events. As in the determination of the time pedestal, these uncertainties 
originate mainly from the random arrival time of cosmic muons relative to the clock cycle. 

The f trig uncertainty of about 10 ns, seen in Fig. [7] (bottom), corresponds to a relative uncer- 
tainty of about 2.5 % on the drift velocity. This is comparable to the fluctuations observed in 
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Figure 8: Mean (top) and width (bottom) Gaussian parameters, as fitted to the distributions of 
the residuals between the reconstructed hits and the reconstructed local segments. The results 
are shown for the r<p super-layers for a representative wheel, after correcting the offset with 
respect to the origin of the residual distribution. Other wheels show similar results. 
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Figure 9: Schematic view of a super-layer section, showing the pattern of semi-cells crossed by 
a track. The quantities t{ in the equations represent the arrival time of the electrons within a 
drift cell. 
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meaning that the residuals calculated with the ^rift anc ^ ^ tri g constants do not represent 



a significant improvement with respect to those shown in Fig.JH] 
The drift velocity distribution measured in one of the VDC chambers (Section 3) during the full 



CRAFT period is shown in Fig. 11 An average drift velocity value of 54.8 pim/ ns is observed, 
which is slightly different (0.5 % higher) from the one obtained from the DT data, mainly be- 
cause of the different shape of the electric field in the DT drift cell. The spread of the distribution 
is better than 0.2 /mi/ ns, and shows that no major variations occurred in the gas mixture or air 
contamination, during the entire data-taking period. 



8 The Calibration Workflow and the Monitoring of the Calibration 
Process 

A fast calibration of the DT system is vital in order to provide the prompt data reconstruction 
with accurate calibration constants. The number of calibration regions is a compromise be- 
tween the need of keeping things simple, not requiring too large event samples, and the need 
of reducing systematic errors by separately calibrating regions where parameters may have 
very different values. As mentioned in previous sections, the super-layer granularity has been 
found to be the most suitable calibration unit. 

In order to reach the precision obtainable with the fast calibration, about 10 4 tracks crossing 
each super-layer are required. During LHC collision and cosmic ray data-taking periods, the 
calibration parameters have to be produced, validated, and made available for use in the re- 
construction within one day of data-taking. However, after the start-up phase, it is anticipated 
that at some point it will no longer be necessary to update the DT calibrations on a daily basis; 
on the other hand, they should be checked against a standard set in order to guarantee their 
stability. 

The workflow of the DT calibration has already been fully embedded into the central CMS cali- 
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Figure 10: Drift velocities computed using the mean-time method for a run with B = T (top) 
and for a run with B = 3.8 T (bottom). Results are shown for the rep super-layers of each 
chamber/ sector of a representative wheel. The other wheels show similar results. 
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Figure 11: Distribution of the drift velocity measured in one Velocity Drift Chamber (VDC) 
during the entire data-taking period. 

bration workflow at a very early stage. A more detailed description of the overall CMS calibra- 
tion and alignment computing workflow used in the CRAFT exercise is given in Refs. fT2llT3l . 
Within the CMS calibration and alignment workflow, particular selections of data, named Al- 
CaReco, were used. They contained a reduced number of events and a reduced event content, 
providing the minimal information to fulfill the requirements of the DT calibration task. The 
sample is saved at the CERN Analysis Facility (CAF) and taken as input to the calibration 
process. The calibration algorithm runs at the CAF and produces a set of constants, which un- 
dergoes a validation procedure before being copied to the central CMS database, where they 
become available to the CMSSW offline software framework. 

The DT calibration workflow has been used also during the Computing, Software, and Anal- 
ysis challenge (CSA08), described in Ref. [fl2|, which simulated with large event samples the 
conditions expected at LHC startup. This exercise simulated the production rate of the calibra- 
tion conditions as it will happen during real collision data-taking. The long CRAFT data-taking 
period served as a thorough test of this workflow with the real detector. 

The quality and stability of the calibration constants is a crucial part of the procedure and 
must be continuously monitored. Therefore, validation procedures have been set up within 
the central CMS Data Quality Monitoring (DQM) framework. A detailed description of the 
CMS DQM structure is given in Ref. |T4|. 

Data quality assessment for the DT calibration constants consists mainly in defining the ac- 
ceptance criteria used to validate the constants, in the monitoring of time stability, and in the 
checking of continuous trends or sudden changes in operating conditions. The quality tests 
to assess the validation of the constants and to monitor their time stability are applied to the 
residual distributions calculated at the different steps of the calibration workflow. The compar- 
ison of the currently produced calibration constants with a reference set gives an indication of 
the stability of each particular calibration constant. 

All the calibration constants described in the present paper have their validation as well as 
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monitoring process, and for each of them detailed and summary DQM plots are provided. The 
DT condition constants have been monitored through the entire CRAFT data-taking period and 
have shown generally a good stability in time. 



9 Drift Velocity Analysis 

The drift velocity obtained with the calibration procedure described in Section 7 is derived from 
the measurements of the drift time and, as already mentioned, is limited by the uncertainty on 
the arrival time of the cosmic ray muons. 
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Figure 12: Computed radial component of the magnetic field in the muon barrel chambers, for 
the different wheels, as a function of z. 



A more detailed analysis of the drift velocity is presented in this section, taking into account 
the precise 3D space-time relationship for the hit reconstruction. In particular, it considers the 
influence of the magnetic field as a function of the position along the wire. 

The presence of a radial magnetic field distorts the drift lines of the drifting electrons, because 



of the Lorentz force, resulting in a variation of the effective drift velocity. Figure 12 shows that 
the radial field component is not very high in the muon barrel chambers, except in the MBl 
chambers of the outer wheels, closest to the endcaps. In these regions, the radial component of 
the magnetic field can be as high as 0.8 T, and changes significantly along the z axis, resulting 
in a variation of the effective drift velocity along the wire of each single cell, for rep super-layers. 
The effect of the magnetic field has been studied in test beams with small prototypes [15J, and 
more recently in the Magnet Test and Cosmic Challenge (MTCC), using cosmic rays in the 
CMS surface hall. These studies showed that the chambers maintain a good trigger and event 
reconstruction functionality, even in the most critical regions [16J. 



In the method presented in this section, a full reconstruction of the trajectory within the muon 
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system is performed to determine the drift velocity. In the first step, a pattern recognition algo- 
rithm is applied to identify hits belonging to the same track. Once the hits have been identified, 
the track is reconstructed under the assumption of a 54.3 jim/ ns nominal drift velocity. In the 
second step the track is refit treating as free parameters the drift velocities at each hit and the 
time of passage of the muon through the chamber. The method is applied to the rep view of 
the track segment in one chamber, where there are eight measured points in most cases. The z 
super-layers, where only four points are available, at most, are less significant for this analysis. 
The drift velocity is taken to be the mean value of the track-by-track drift-velocity distribution. 
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Figure 13: Mean values of the drift velocity for the MB2 chambers of Wheel +2, using the 
mean-time (squares) and fit (circles) methods. The differences between sectors when using the 
mean-time method are due to £ t rig uncertainties that are not present in the fit method. 



Figure 13 shows the mean values of the drift velocity for the MB2 chambers of Wheel +2, using 
the mean-time method described in Section 7 and the fit method described here. When using 
the mean-time method, the drift velocities have large systematic fluctuations from one sector 
to the other. This is related to the errors on the £ t rig determination described in Section 6, which 
cancel when the fit method is used. 



The average drift-velocity values from the fit method, for all the chambers, are shown in Fig. 14 
for runs without and with magnetic field (of 3.8 T). 

The data at B = T show an average value of 54.5 /mi/ ns for the drift velocity and a standard 
deviation indicating that differences between chambers are in the order of 0.2 %. For B = 3.8 T, 
a second peak is observed at 53.6 f/m/ ns. This peak corresponds to the MB1 chambers of the 
external wheels (Wheel +2 and Wheel —2) and is due to the presence of a higher radial magnetic 
field. 



Similar values of the drift velocity have been obtained using the same calibration procedure 
applied to the simulated pp collision data. These results, presented in Ref. [17J, indicate that 
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Figure 14: Drift velocities for B = T (left) and B = 3.8 T (right). The small peak on the right 
panel corresponds to the MB1 chambers of Wheel +2 and Wheel —2, and shows the influence 
of a higher magnetic field in these regions. 

the calibration algorithm delivers a more uniform response in the case of collision data and that 
a large fraction of the fluctuations observed in the drift velocity calibration from CRAFT data 
may be attributed to the topology and timing of the cosmic ray events. 
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Figure 15: Drift velocities calculated using the muon track fit method described in this section. 
The values are shown as a function of the z position (measured by the z super-layers), and for 
B = T and B = 3.8 T. 



The effect on the drift velocity of the variation of the radial magnetic field along the z coordinate 
is shown in Fig.[l5j as calculated with the fit method. Positive wheels (+1 and +2) are not in the 
figure but show the same behavior as their symmetric wheels (—1 and —2, respectively). 

The presence of the radial component of the magnetic field affects, as expected, only the MB1 
chambers, primarily in the external wheels but some effects are also observed in Wheels +1 and 
— 1. The variation along z for the MB1 chambers of Wheels +2 and —2 is below 3 %, less than 
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expected from the MTCC results [16], after taking into account the differences of the magnetic 
field conditions between both periods (B = 4 T during the MTCC in the surface hall, B = 3.8 T 
in the underground experimental hall). 

This analysis of the drift velocity is very sensitive to the field strength and, in fact, provided 
the first evidence of a systematic deviation from the true field strength in the field map. A 
new detailed calculation of the magnetic field has been performed [TlBI , and the new analyses, 
currently in progress, show a z dependence in satisfactory agreement with the expectations. 



All sectors in the same wheel show the same behavior, as illustrated in Fig. 16 where the values 
of the drift velocity along the z axis are shown for the MB1 chambers of some representative 
sectors of Wheels +2 and 0. 

The drift velocity calculation, performed in this section, provides a better spatial resolution of 
the chambers with respect to the one obtained in Section 7. This improvement is obtained with 
an extended track fit method which determines the drift velocity and the time of passage of the 
muon simultaneously with the regular track parameters. The detailed analysis of the spatial 
resolution for the cosmic ray data taking in 2008 is given in Ref. |19| . The value obtained is 
about 250 /mi, in fair agreement with the requirements for collision dataa [3] . 



10 Summary 

This paper describes the calibration of the CMS Drift-Tubes system and presents results from 
the cosmic ray data-taking period which took place in 2008. 

The complete calibration workflow has been applied to the data. It performed efficiently, mon- 
itoring the stability of the produced constants, and delivering with very low latency the cali- 
bration constants to the conditions database used by the offline reconstruction. 

The first calibration step is the identification and masking of noisy channels to have a clean 
structure of the drift time distribution. The fraction of noisy cells was stable and about 0.01 %. 
The average noise rate was ~ 4 Hz. 

The time pedestals, after having been corrected for the inter-channel synchronization, noisy 
channels, and the time of flight between upper and lower sectors, show a constant behavior 
in the entire DT system. Due to the particular topology of the cosmic ray events, the time 
pedestals are poorly defined for the sectors with chambers in the vertical plane, where cosmic 
ray tracks with large impact angles are measured. For all the other sectors, an uncertainty of 
the order of 10 ns is observed. This value agrees with the uncertainty of the arrival time of 
cosmic ray muons within the clock cycle. 

The drift velocity calibration results show an approximately constant value of 54.3 jim./ ns for 
all the chambers of the DT system, with a relative systematic uncertainty of 2.5 %. This uncer- 
tainty originates from the measured drift time, used in the mean-time method, which is limited 
by the uncertainty of the arrival time of cosmic ray muons. This explains why the obtained 
spatial resolution is worse than would be expected with collision data. 

A more refined analysis of the drift velocity has been performed, exploiting the full potential 
of the CMS offline software for data reconstruction. It uses a track fitting procedure which 
leaves as free parameters the drift velocity and the time of passage of the muons through the 
chambers. Cosmic ray data with and without magnetic field have been studied. Without mag- 
netic field, a constant average value of 54.5 }im/ ns has been observed, with an error of 0.2 %; 
when the field strength is 3.8 T, the innermost chambers of the external barrel wheels measure 
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Figure 16: Drift velocities as a function of the local z position for MB1 chambers of some rep- 
resentative sectors of Wheel (top) and Wheel +2 (bottom). Different sectors are indicated by 
different grey tones. 
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a lower value, as expected, of about 53.6 Jim/ ns. These results confirm what was observed in 
an analysis performed on simulated collision data and provide a spatial resolution that is close 
to the design performance. 
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